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1 Introduction 



The physics of high-density QCD has become an increasingly active subject 
of research, both from experimental and theoretical points of view. Presently, 
and in the near future, the collider facilities such as the DESY ep collider 
HERA {ep, eA), Fermilab Tevatron {pp, pA), BNL Relativistic Heavy Ion 
Colhder (RHIC), and CERN Large Hadron Collider (LHC) {pp, AA) will be 
able to probe new regimes of dense quark matter at very small Bjorken x 
or/and at large A, with rather different dynamical properties. The descrip- 
tion of these processes is directly associated with a correct description of the 
dynamics in this kinematical region. 

The behavior of the cross sections in the high energy limit {s oo) and 
fixed momentum transfer is expected to be described by the BFKL equation 
[0]. The simplest process where this equation applies is the high energy 
scattering between two heavy quark- ant iquark states, i.e. the onium-onium 
scattering. For a sufficiently heavy onium state, high energy scattering is a 
perturbative process since the onium radius gives the essential scale at which 
the running coupling as is evaluated. This process was studied in the dipole 
picture p|, ||, §], where the heavy quark- antiquark pair and the soft gluons 
in the limit of large number of colors Nc are viewed as a collection of color 
dipoles. In this case, the cross section can be understood as a product of 
the number of dipoles in one onium state times the number of dipoles in the 
other onium state times the basic cross section for dipole-dipole scattering 
due to two-gluon exchange. The cross section grows rapidly with the energy 
{a oc a^gC^p-i)^, where {ap - 1) = ^^^ln2 and Y = Ins/Q^) because 
the number of dipoles in the light cone wavefunction grows rapidly with the 
energy. In Mueller demonstrated that the QCD dipole picture reproduces 
the BFKL physics. 

One of the main characteristics of the BFKL equation is that it predicts 
very high density of partons in the small x region. Therefore a new dy- 
namical effect associated with the unitarity corrections is expected to stop 
further growth of the parton densities. The understanding of the unitarity 
corrections has been a challenge of perturbative QCD (PQCD). 

About seventeen years ago, Gribov, Levin, and Ryskin (GLR) 0] per- 
formed a detailed study of this problem in the double logarithmic approxi- 
mation (DLA). They argued that the physical processes of interaction and 
recombination of partons become important in the parton cascade at a large 
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value of the parton density, and that these shadowing corrections could be 
expressed in a new evolution equation - the GLR equation. This equation 
considers the leading non-ladder contributions: the multi-ladder diagrams, 
denoted as fan diagrams. The main characteristics of this equation are that 
it predicts a saturation of the gluon distribution at very small x, it predicts a 
critical line, separating the perturbative regime from the saturation regime, 
and it is only valid in the border of this critical line. Therefore, the GLR 
equation predicts the limit of its validity. In the last decade, the solution 
[^, 13, ^ and possible generalizations |lOl of the GLR equation have been 
studied in great detail. Recently, an eikonal approach to the shadowing cor- 
rections was proposed in the hterature |jTl], |12|, [l^. The starting point of 
these papers is the proof of the Glauber formula in QCD which con- 
siders only the interaction of the fastest partons with the target. In a 
generalized equation which takes into account the interaction of all partons 
in a parton cascade with the target in the DLA limit was proposed by Ay- 
ala. Gay Ducati, and Levin (AGL). The main properties of the generalized 
equation are that (i) the iterations of this equation coincide with the iter- 
ation of the Glauber-Mueller formula; (ii) its solution matches the solution 
of the DGLAP evolution equation in the DLA limit of PQCD; (iii) it has 
the GLR equation as a limit, and (iv) contains the Glauber-Mueller formula. 
Therefore, the AGL equation is valid in a large kinematic region. 

The AGL equation resums all multiple pomeron exchanges in the DLA 
limit. Its assintotic solution is given by xG oc Q"^ In {1 / x) , where R is the 
size of the target, i.e. differently from the GLR equation, it does not predict 
saturation of the gluon distribution in the very small x limit. 

The unitarity corrections in the leading logarithmic limit can be studied 
using the dipole picture. In this picture the unitarity corrections (multi- 
ple interactions between the onia) become important at high energies [|T5| , 
P^ . These corrections in general are neglected as being of higher order in 
1/N(.. However, when the energy is high, there is a large number of dipoles 
in each onium, and the total number of possible interactions is of order 
g(ap-i)(j/+(y-j/)) i^^Yie product of the number of dipoles in each onium) 
So when a^e*^"^"^-'^ ~ 1, one should take into account the multiple scat- 
tering corrections - multiple pomeron exchange - despite the fact that they 
are suppressed by l/N^. Recently, an equation which includes all pomeron 
exchanges in the leading logarithmic approximation using the dipole picture 
was proposed (Eq. (15) in Ref. [O])- In this paper we will denote this 
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equation as K equation. As the K equation and the AGL equation resums 
the same group of diagrams, both should be identical in a common limit. 
Our goal in this paper is to demonstrate that the K equation reproduces the 
AGL equation and, consequently, the GLR equation in the DLA limit. 

This paper is organized as follows. In Section 2 the AGL approach for the 
unitarity corrections is briefly reviewed. We clarify some steps and obtain 
the GLR equation as a limit case. In Section 3 we obtain the nuclear struc- 
ture function F2 in the Glauber-Mueller approach and compare it with the 
expression proposed in . It allows to obtain precisely the relation between 
the propagator of the qq through the nucleus and the gluon distribution. In 
Section 4, we obtain the AGL equation from the dipole picture and, as a 
limit, the GLR equation. Finally, in Section 5, we present our conclusions. 



2 The AGL Equation 

In the nucleus rest frame we can consider the interaction between a virtual 
colorless hard probe and the nucleus via a gluon pair [gg) component of the 
virtual probe. In the region where x « l/2mR [R is the size of the target), 
the gg pair crosses the target with fixed transverse distance rt between the 
gluons. Moreover, at high energies the lifetime of the gg pair may substan- 
tially exceed the nuclear radius. The cross section for this process is written 




as 



JO J IT 

where G* is the virtual colorless hard probe with virtuality Q^, z is the 
fraction of energy carried by the gluon and is the wave function of the 
transverse polarized gluon in the virtual probe. Furthermore, aa9+^{z,rf) is 
the cross section of the interaction of the gg pair with the nucleus. 

As our goal in this paper is to obtain the AGL from the dipole picture, we 
should make a transformation of the above equation (|l]) for qq dipoles, since 
in this picture they are the basic configuration. Considering that a^^"'""^ = 
{CA/CF)a'^^ we have that the cross section for the interaction between the 
virtual probe G* and the nucleus via the qq component of the virtual probe 
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is given by 



G*A 



Ca 
Cf 



dz 



2nq+Ai 



(2) 



To estimate the unitarity corrections we have to take into account the 
rescatterings of the quark- ant iquark pair inside the nucleus. The contribu- 
tions of the rescatterings can be estimated using the Glauber-Mueller ap- 
proach proposed in ref. [|n|]. Following the same steps used in |TT| , i. e. 
considering the s-channel unitarity and the eikonal model, one obtains that 
the a^*"^ cross section is written as 



G*A 



Cp Jo 



TT 



^fiQ^\n,z)\■'2[l-e 



(3) 



where x' = x/ {z rf Q"^) {x is the Bjorken variable), bt is the impact parameter, 

S{bt) = {A/7cR'^)e is the gaussian profile function and is the cross 
section of the interaction of the qq pair with the nucleons inside the nucleus. 

^{3a^)/A)7r'r^xG^ix,^), 



In HTSi the authors have shown that cr^ 



where xGn{x, t) is the nucleon gluon distribution. 

The wavefunction was calculated in |lT] using the technique of 
Ref. |T^. Here we only explicitate the result for the squared wavefunction, 
which is given by 



^f{Q'',n,z) 



zl 



n 



(4) 



where = Q'^z{l — z) and the Ki are the modified Bessel functions. The 
main contribution in expression comes from the region of small z, where 
ert <C 1. Using the expansion of the modified Bessel functions for small 
values of the argument, results that the squared wavefunction simplifies to 



|vl>^ [Q\rt,x,z)\' = — 

zrf 



The condition ert <^ 1 implies that 



(5) 



(6) 
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where the last inequahty comes from symmetric pairs {z = 1/2). From the 
expression (j^) results that x' > x. 

The relation a'^'^{x,Q^) = {4:'iT^as/Q'^)xGA{x,Q'^) is valid for a virtual 
probe G* with virtuahty Q^. Consequently, using the expression (H) of the 
squared wavefunction and making the change of variables z ^ x' , we obtain 
that the Glauber-Mueller formula for the interaction of the qq pair with the 
nucleus is written as 

xGa{x,Q ) = —— — I —J ^2[l-e ]. (7) 



n"^ Cp Jx x' J-4 Tcrf J vr 



The lower limit in integration comes from the expression (P). 

The space-time picture of the Glauber-Mueller approach is presented in 
Fig. ||. It takes into account only the interaction of the fastest partons with 
the target. As in QCD we expect that all partons from the cascade interact 



with the target, a generalized equation was proposed in |jTT . 

The AGL equation can be obtained directly from the Eq.(|^ differentiat- 
ing this formula with respect to Inl/x and /nQ^/Ag^^. Therefore the AGL 
for the interaction of qq dipole is given by 



d'xGAix.Q') 2Q'Ca fd% ,^ _i 



dln{l/x)dln{QyAlcD) Cp J n 



[1 _ g-ta^'(x,Q2)5(feO] ^ (8) 



where the dependence of cr^ in the virtuality of the virtual probe results from 
the derivative. The nonperturbative effects coming from the large distances 
are absorbed in the boundary and initial conditions. This equation is valid 
in the double logarithmic approximation (DLA). 

For a central collision (6 = 0) the AGL equation (||) reduces to 



d^xGAix,Q^) Ca2Q^R^ 



dln{l/x)dln{Q^/KlcD) Cp 



TT 



2 



[1 _ ^-hVM-)Sio),^ ^ (9) 



where we have set dh"^ = for a 6 = collision. This limit is important since 



the cross section is strongly unitarized at small impact parameters |T^, |T6 

Considering that the transverse cross-sectional area of the nucleus is S± = 
nR^ and that S'(O) = A/i^nR"^), the AGL equation for 6 = is obtained as 

. = - e-^*-^'-«''] . (10) 

dln{l/ x)dln{Q^ / Kf^CD) ^r^ 
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where Nc = 3 and we have assumed that Cp = Nc/'i in the large Nc hmit. 



The Eq. (|T^) agrees with the expression (2) of . 

The AGL equation takes into account the interaction of all partons in 
a parton cascade with the target. In other words, the AGL equation takes 
into account that each parton in the parton cascade interacts with several 
nucleons within of the nucleus (Glauber multiple scattering). In Fig. |^ we 
present the result of the first iteration from the Glauber-Mueller formula, 
which is one of the terms summed by the AGL equation. 

The GLR equation can be obtained directly from the Eq. ([To|). If we 
expand the right hand side of this equation to the second order in xGa we 
obtain 



dln{l/x)dln{QyAlcD) ^ ^v'^/ q 

which is the GLR equation for a cylindrical nucleus case (See Eq. (19) in 
[p^ ). Moreover, if the unitarity corrections are small, only the first order in 
xGa contributes. In this limit the equations ([lOD and (|lTD matches with the 
DGLAP evolution equation in the DLA limit. 

Therefore the AGL equation (i) matches the DGLAP evolution equation 
in the DLA limit of PQCD; (ii) it has the GLR equation as a limit, and (iii) 
contains the Glauber-Mueller formula. We have that the AGL equation is 
valid in a large kinematic region. 

Recently, a comprehensive phenomeno logical analysis of the behavior of 
distinct observables was made for the HERA kinematical region using the 
Glauber-Mueller approach p!3| , |2^ , p^ . In this kinematical region the solu- 



tions from the AGL equation and the Glauber-Mueller approximately coin- 
cide. The results from these analysis agree with the recent HERA data and 



allows to make some predictions which will be tested in a near future. In 
we have analysed the behavior of the longitudinal structure function Fl and 
the charm component of the proton structure function and have shown 
that our results agree with the HI data. New data, with better statistics. 



will allow to evidentiate the unitarity corrections. In p3| we have shown that 
the recent ZEUS data can be described if the unitarity corrections for the F2 
slope are considered. Our main conclusion is that the unitarity corrections 
cannot be disregarded in the HERA kinematical region. 
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3 The Nuclear Structure Function 



The unitarity corrections to the nuclear structure function can be estimated 
in the rest frame of the target. This intuitive point of view was proposed 



by V. N. Gribov many years ago pj]. Gribov's assumption is that at small 
values of x the virtual photon fluctuates into a qq pair well before the in- 
teraction with the target, and this system interacts with the target. This 
formalism has been established as an useful tool for calculating deep inelas- 
tic and related diffractive cross section for 7*p scattering in the last years 



2^, The Gribov factorization follows from the fact that the lifetime of 
the qq fluctuation is much larger than the time of the partonic interactions. 
According to the uncertainty principle, the fluctuation time is ~ — , where 
m denotes the target mass. 

The space-time picture of the DIS in the target rest frame can be viewed 
as the decay of the virtual photon at high energy (small x) into a quark- 
antiquark pair long before the interaction with the target. The qq pair sub- 
sequently interacts with the target. In the small x region, where x -C 
the qq pair crosses the target with fixed transverse distance rt between the 
quarks. It allows to factorize the total cross section between the wave func- 
tion of the photon and the interaction cross section of the quark- ant iquark 
pair with the target. The photon wave function is calculable and the inter- 
action cross section is modelled. Therefore, the nuclear structure function is 
given by 

F,^{x,Q') = JdzJ ^\^{z,r,)\'a^-^^\z,n) , (12) 

where 

l^(^,r.)P = ^E4{[-' + (l-^)Vi^i(^r.r + ^?i^o^(er.)n , (13) 

ttem is the electromagnetic coupling constant, e = z{l — z)Q'^ + rrij, m/ is 
the quark mass, n/ is the number of active flavors, is the square of the 
parton charge (in units of e), i^o,i are the modified Bessel functions and z 
is the fraction of the photon's light-cone momentum carried by one of the 
quarks of the pair. In the leading log(l/a:) approximation we can neglect the 
change of z during the interaction and describe the cross section a'^^~^^{z, rf) 
function of the variable x. 
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We estimated the unitarity corrections considering the Glauber multiple 
scattering theory |27], which was probed in QCD [14]. The nuclear collision 
is analysed as a sucession of collisions of the probe with individual nucle- 
ons within the nucleus, which implies that the F2 structure function can be 



written as 11 



J J Ti J 71 



} (14) 



where the opacity flqq{x, rt, ht) describes the interaction of the qq pair with 
the target. 

In the region where Vtqq is small {^qq ^ 1) the ht dependence can be 



factorized as Vt 



VtqqS{ht 



with the normalization / cPht S{ht) = 1. 



The eikonal approach assumes that the factorization of the bt dependence 
flqq = flqqS{bt), whlch Is Valid in the region where Qqq is small, occurs in the 
whole kinematical region |jl2| . The main assumption of the eikonal approach 
in pQCD is the identification of opacity flqq with the gluon distribution. In 
[|T^ the opacity is given by ilqq = afj. 

For a central collision {h = 0) the Eq. (|l^ reduces to 



-W / dz 



d^n 



*(z,rt)|" 2{l-e-^"^'^'^"^} . (15) 



5(0)1 



Substituting cr^ and 5(0) we obtain that the nuclear structure function for 
6 = is given by 



-W / dz 



d^n 



7T 



^{z,rt)\^ 2{l-e 



-r^AxG{x,l/r'^) 



} m 



This expression estimates the unitarity corrections to the nuclear struc- 
ture function for central collisions {b = 0) in the DLA limit using the Glauber- 
Mueller approach. 

Comparing the Eq. ([T6|) with the expression to the nuclear structure 
function in the dipole picture proposed in [|l7l, we obtain the representa- 
tion for the total cross section of the qq pair interacting with the nucleus, 
A^(a;oi, bo = 0,Y), used in that reference. 



iV(xKi,6o = 0,r)= 2{l-e 



} 



(17) 
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where xqi = xq — xi = rt [xq (xi) is the transverse coordinate of the 
quark (antiquark)] and Y = ln{s/Q'^) = ln(l/x). This relation between 
N{x~Qi, bo = 0, Y) and the gluon distribution is vahd in the DLA hmit and is 



two times the propagator of the qq pair through the nucleus obtained in ^0 
at large (See Eq. (6b) of |T 



We will assume that the correct expression for N in the DLA limit is 
given by (|T^. In the next section we will use this expression as an input in 
the evolution equation for N obtained in the dipole picture. 



The AGL Equation from the Dipole Pic- 
ture 



Considering the multiple pomeron exchange, Kovchegov has obtained an 
evolution equation for A^(a;oi, bo, Y) in the leading logarithmic approximation 



N{xoi,bo,Y) = -7(xoi,6o)exp 



-^In f ^ ] Y 



TC 



P 



TT^ JO 



dy exp 



^ln(j-\ (Y-y) 



X / (fx2 2 °\ [2A^(a;o2, &o, v) - N{xo2, &o, y)N{xi2, bo, y)] 



(18) 



where Xij = Xi — xj is the size of the dipole with a quark in the transverse 
coordinate Xi and the antiquark in Xj, 7(a;oi,6o) is the propagator of the qq 
pair through the nucleus and p is an ultraviolet cutoff in the equation which 
disappears in the physical quantities. We have reobtained this equation and 
denote it as the K equation. 

The equation (|1^) was obtained considering the scattering of a virtual 
photon with a nucleus. The physical picture for this interaction is the same 
as the Glauber-Mueller approach. The incoming virtual photon generates 
a qq pair which develops a cascade of gluons, which then scatters on the 
nucleus. In the large Nc limit the gluon can be represented as a qq pair. 
Therefore, in this limit and in the leading logarithmic approximation, the 
cascade of gluons can be interpreted as a dipole cascade, where each dipole in 
the cascade interacts with several nucleons within the nucleus. Therefore, as 
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the K equation and the AGL equation, although with distinct basic objects, 
resums the multiple scatterings of its respective degree of freedom, we expect 
that both coincide in a common limit. 

In the double logarithmic limit, where the momentum scale of the photon 
is larger than the momentum scale of the nucleus Kqcd-, we take the large 
limit from ([T^), which reduces to 



dY ~ ^;~^oiL 7^-^2 



A. I. f^[2N{xo2,bo,Y) 

02) 

N{x-o2,bo,Y)N{xo2,bo,Y)] . (19) 



This equation considers the evolution of the dipole transverse sizes from a 
small scale Xqi to a large scale 1/Aqcd- Differentiating the above expression 
with respect to In (l/a^oiAQj-^^,) a double differential equation for A^(xoi, bo, Y) 
follows 

[2 - N{xoi,bo, Y)\N{xoi, bo, Y) . (20) 



dYdln{l/xl,AlcD) 



TT 



The physical picture for the dipole evolution in the DLA limit is that the 
produced dipoles at each step of the evolution have a much greater transverse 
dimensions that the parent dipoles . 

In the connection between the quantity N{xoi, bo, Y) and the nuclear 



gluon distribution was discussed. As there is some freedom in the definition 
of the gluon distribution, the choice for the connection between the two 
functions was arbitrary. Here we use the result obtained in the previous 
section [Eq. (17)] to make more explicit this connection. If the Eq. (|T^) is 



expanded to the first order in xGa, our result differs from the choice used in 
[ p!7[] by a factor of 2. 

Substituting Eq. (|I^ into Eq. (pO|) and using xqi ^ 2/Q |T^, which is 
valid in the double logarithmic limit, we end up with 

d'xGAix,Q') _ N,CpS^_^2^^_^.^^^.GA^,Q^)^ ^ ^21) 



dln{l/x)dln{QyAlcD) 



which exactly corresponds to the AGL equation [Eq. (|TOD]- As demonstrated 
in the Section ^ the GLR equation is a straithforward consequence from the 
AGL equation. 
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A last comment is important. Using the relation between N and xG from 

[p!7[l a factor e remains in the right-hand side from the Eq. (|2T|) . 

Therefore, the AGL equation only could be justified for a configuration of 



very small dipoles, as was discussed in [|17|. Consequently, the connection 
between the dipole and AGL approaches deals on the correct physical relation 
between the function and the nuclear gluon distribution, given here by Eq. 

a- 



5 Conclusions 

The AGL equation, obtained using the Glauber-Mueller approach, resums 
all multiple pomeron exchanges in the double logarithmic limit. Its solution 
allows to make predictions which agree with the recent HERA data. In 
this paper we have analysed another equation which includes all multiple 
pomeron exchanges proposed by Kovchegov recently. The K equation was 
obtained in the leading logarithmic approximation using the dipole picture. 
We agree with this result. However, as the relation between the function 
A^, which evolution is described by the K equation, and the nuclear gluon 
distribution has some freedom, we obtained this connection considering the 
nuclear structure function in the Glauber-Mueller approach. Substituting 
this relation in the evolution equation for A^ we have shown that it reduces 
directly to the AGL equation and, as a limit, to the GLR equation in the 
double logarithmic limit. This result shows that the AGL equation is a good 
candidate for the unitarized evolution equation at small x in the DLA limit, 
supported by two different frameworks describing small x phenomena. 

Another candidate for the unitarized evolution equation was proposed by 



Jalilian-Marian et al. [|2^. These authors have derived a general evolution 
equation for the gluon distribution in the limit of large parton densities and 
leading logarithmic approximation, considering a very large nucleus. This 



work is based on effective Lagrangian formalism for the low x DIS and 
the Wilson renormalization group. In the general case the evolution equation 
is a very complicated equation, which does not allow to obtain analytical 
solutions. Recently, these authors have considered the DLA limit on their 
result pO| and have shown that the evolution equation reduces to an equation 



with a functional form similar, but no identical, to the AGL equation. We 
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believe that a more detailed analysis of the approximations used in both 
equations will allow to demonstrate the equivalence of both equations in a 
common limit BT . 
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Figure 1: Space-time picture of the Glauber-Mueller approach. The virtual 
colorless hard probe G* decays into a pair with transverse distance which 
interacts with the n nucleons within the nucleus. 



Figure 2: First iteration of the Glauber-Mueller approach. The AGL equa- 
tion takes into account the interation of all partons of the parton cascade 
with several nucleons within of the nucleus. 
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